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Abstract. The half lives of proton radioactivity of proton emitters are investigated theoretically. Proton- 
nucleus interaction potentials are obtained by folding the densities of the daughter nuclei with a finite 
range effective nucleon-nucleon interaction having Yukawa form. The Wood-Saxon density distributions for 
the nuclei used in calculating the nuclear as well as the Coulomb interaction potentials are predictions of 
the interaction. The quantum mechanical tunneling probability is calculated within the WKB framework. 
These calculations provide reasonable estimates for the observed proton radioactivity lifetimes. The effects 
of neutron-proton effective mass splitting in neutron rich asymmetric matter as well as the nuclear matter 
incompressibility on the decay probability are investigated. 
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1 Introduction 

In recent years new half life measurements have been per- 
formed in order to have a better understanding of the pro- 
ton and alpha decay processes in the region of proton-rich 
nuclei pQ. For these nuclei the Q value for proton emis- 
sions is positive and therefore there is a natural tendency 
to shed off excess protons. These data are very useful for 
the analysis of possible irregularities in the structure of 
these proton- rich nuclei 0[3]. They are also of great in- 
terest in rapid proton capture nucleosynthesis processes. 
Some new results for proton radioactivity in this region of 
proton-rich nuclei have indicated that the proton emission 
mode is rather competitive with the alpha decay process 
[3UH[S] . Proton radioactivity may be used as a tool to ob- 
tain spectroscopic information because the decaying pro- 
ton is unpaired in the orbit. These decay rates are sensitive 
to the Q values and the orbital angular momenta which 
in turn help to determine the orbital angular momenta of 
the emitted protons. 

Since the observation of proton radioactivity is com- 
paratively recent, several theoretical approaches that have 
been employed to study this exotic process, such as the 
distorted- wave Born approximation 6 , the density depen- 
dent M3Y (DDM3Y) effective interaction the effec- 
tive interaction of Jeukenne, Lejeune and Mahaux (JLM) 
[8], the unified fission model [3], the coupled-channels ap- 
proach [TU] and the effective/generalized liquid drop mod- 
els [HlHl?] are also quite recent [TBI 14) . In the present 
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work, quantum mechanical tunneling probability is cal- 
culated within the WKB framework using proton-nucleus 
interaction potentials obtained from folding the density of 
the residual daughter nucleus with a finite range effective 
nucleon-nucleon interaction having a single Yukawa term 
(YENI) in the finite range part [TSlfTB] . These calculations 
provide good estimates for the observed proton radioactiv- 
ity lifetimes. In the present calculation we shall examine 
the effects pf neutron-proton (n-p) effective mass splitting 
as well as the nuclear matter incompressibility, K(flo), on 
the decay probability of the proton emitters. The n-p effec- 
tive mass splitting is connected to the momentum depen- 
dence aspect of neutron and proton mean fields in asym- 
metric nuclear matter (ANM). Theoretical predictions of 
different models on this important issue can be divided 
into two distinct groups depending on whether the neutron 
effective mass goes above the proton one [17,18,19,2(31 1211 
[2"2"] or the other way around [2ifll241l25l[2"o1[2"71 . Experimen- 
tal as well as theoretical attempts [281129] to resolve the 
problem has not been successful as yet. In the work in 
Rcf.JTB] it is shown that if the finite range exchange in- 
teraction acting between an unlike nucleon pair, v™,{f), 
is stronger compared to the finite range exchange interac- 
tion, v ex , between a pair of like nucleon, the neutron ef- 
fective mass will be predicted to go over the proton one in 
neutron rich ANM. On the other hand if v l ex (r) is stronger 
compared to v™ x (r) then proton effective mass will go over 
the neutron one. The results of Dirac-Brueckner-Hartrcc- 
Fock (DBHF) calculations jTM2"T] along with the decreas- 
ing trend of Lane potential extracted from experimental 
data on nucleon-nucleus scattering and reaction [3T3113T1 
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32 . 33] has led to accept amongst a larger community that 
neutron effective mass goes above the proton one in neu- 
tron rich matter although the controversy is not yet com- 
pletely resolved. Under the circumstance the magnitude 
of effective mass splitting remains as an open problem as 
different models give widely divergent results. Similarly, 
the nuclear matter incompressibility is another important 
quantity whose value ranges between 240 ± 20 MeV as has 
been estimated from studies of isoscalar giant monopole 
resonances in nuclei. In the present work we shall also ex- 
amine the possible effect of the variations of these nuclear 
matter parameters on the calculated proton half-lives of 
the proton emitters. 

In section 2 the calculation of proton-nucleus (p-N) 
interaction potential for any general effective interaction 
has been discussed. The formalism has been extended to 
the calculation with the YENI. The determination of the 
parameters required in the study of the nuclear matter is 
briefly discussed. The fixation of the free parameter of the 
interaction along with the Wood-Saxon (WS) density dis- 
tribution for the nuclei have been obtained by adopting 
a simultaneous minimization procedure to reproduce the 
binding energies of nuclei. Contributions of the different 
parts of YENI to the p-N nuclear potential and the self- 
consistent evaluation of the finite range exchange part are 
provided. In section 3 WKB tunneling procedure for cal- 
culation of decay probability of emitted proton has been 
discussed. The last section contains discussions of the re- 
sults obtained and conclusions. 



k{R) = sJ-^{E cm - V N (R) - V C {R)) (2) 

where E cm , Vn(R) and V C (R) are the center of mass en- 
ergy, p-N interaction potential and Coulomb potential at 
the center of mass R, respectively. The center of mass 
and relative coordinates are given by R = (r + r')/2 and 
t = (r — r') respectively. It may be seen from Eq.(l) 
that in the calculation of Vn(t) the knowledge of k(R) 
is required in which Vjv(i?) appears and hence requires 
a self-consistent calculation. The total interaction energy 
between the proton and the residual daughter nucleus 
V(r) = V N (r) + V c (r) + h 2 l(l + l)/(2pr 2 ), the sum of 
the nuclear interaction energy, the Coulomb interaction 
energy and the centrifugal barrier where I is the angular 
momentum carried away by the proton-daughter nucleus 
system. Here p = M p Md/MA is the reduced mass, M p , 
Md and Ma are the masses of the proton, the daughter 
nucleus and the parent nucleus respectively, all measured 
in the units of MeV/c 2 . 



2.2 Simple finite range effective interaction and the 
proton-nucleus potential 

The simple parameterization of finite range effective in- 
teraction |16j used in this work for calculating proton ra- 
dioactivity of the spontaneous proton emitters is given by, 



2 The proton-nucleus interaction potentials 

2.1 The folded proton-nucleus interaction potential 

The proton-nucleus potential is obtained by folding the 
density distribution of the nucleus over the interaction of 
the incident proton with nucleons of the nucleus. It is given 

by M 

V N (r) = J[ Pp (r')v p /(\r - r'\) + Pn (r f )v^(\r - r'|)]dV 
+ Jp p (ry)j (k(R)\r-r'\)v^(\r-r'\)d 3 r' 
+ J p n (r,r')j (k(R)\r - r'\)vg(\r - r'|)dV 
+ rearrangement terms (1) 

where r and r' are the distances of the incident proton 
and the nucleon of the nucleus, respectively, from the ori- 
gin taken at the center of the nucleus. The last term in 
Eq.(l) is the rearrangement term that arises from the 
explicit density dependence of the effective interaction. 
Pi(r,r'),i — p,n are the density matrices that take non- 
local effects into account, and v d / ex is the direct/exchange 
part of the effective interaction averaged over space, spin 
and isospin of both the interacting nucleons. jo is the ze- 
roth order spherical Bessel function. k(R) is the wave 
number of the incident proton at the center of mass R 
of the incident proton and nucleon of the nucleus and is 
given by, 



v eff (r - r') = *q(1 + x P a )5(r - r') 

+ (W + BP a - HP T - MP a P T )f a {\r - r'\) (3) 

where / Q (|r — r'Q, is a short range interaction of conven- 
tional form, such as, Yukawa, Gaussian or exponential and 
specified by a single parameter a, the range of interaction. 
This effective interaction contains altogether eleven ad- 
justable parameters, namely, to, xq, ts, x%, b, 7, W, B, H, M 
and a. P a = (l + £r 1 .o- 2 )/2 and P T =(1 + T\.t%)/2 are the 
spin and isospin exchange operators respectively. This in- 
teraction has been used in the studies of momentum and 
density dependence of both symmetric and asymmetric 
nuclear matter at zero and finite temperatures [15,16 
as well as in the calculation of bulk properties of neu- 
tron stars [35] and equation of state (EOS) of beta stable 
n+p+e+p matter, i.e., neutron star matter(NSM) [55]. In 
these studies we require a total of nine parameter combina- 
tions, namely, a, 6, 7, e l , £q', e~ l , e l ex and e^ x out of the 
total eleven interaction parameters for the complete de- 
scription of asymmetric nuclear matter and their relations 
to the interaction parameters are given in Ref. [35 . Out 
of these nine parameters required for a complete descrip- 
tion of ANM only six, namely a, b, 7, (e l + Eq 1 ), (e l + s" 1 ) 

and (e l ex + e™ x ) are required to describe the EOS of sym- 
metric nuclear matter (SNM). The careful adjustment of 
these six parameters so as to provide a correct momentum 
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dependence of the mean field as well as density depen- 
dence of the EOS in SNM is discussed in the Refs. [35ll37] . 
The crucial advantage of the procedure adopted to con- 
strain these six parameters in SNM is that the momen- 
tum dependence of the mean held can be varied with out 
changing the density dependence of the EOS of SNM and 
vice- versa is also true. The momentum dependence of the 
mean held in SNM is decided by the finite range exchange 
strength parameter (e l ex + e"£) and the range a, whereas, 
the stiffness of the EOS is determined by the parameter 7 
in the exponent. Under the consideration that the inter- 
action between pairs of like (I) and unlike (ul) nucleons 
have same range but differ in strength, the study of ANM 
now, requires the correct splittings of the three parame- 
ters (£q + £q 1 ), (e' + e" z ) and (e l ex + e^ x ) into two specific 
channels for interaction between pairs of like and unlike 
nucleons. In absence of adequate constraints, either exper- 
imental or theoretical, to decide the splitting of these three 
strength parameters the procedure that we have adopted 
in our study of ANM has been discussed in Refs. [161155] . 
The splitting of finite range exchange strength parameter 
(s l ex + into like and unlike channels decides the n- 
p effective mass splitting in ANM. The possible range of 
splitting into the like channel, i.e., e ex> can be from to 
(e l ex + and accordingly the exchange strength in the 
unlike channel e"f x is decided. For e l ex in between and 
(e l ex + e^ x )/2 the neutron effective mass is predicted to lie 
over the proton one and for e l ex in the range {e l ex + s^ x )/2 
and (e l ex + e^ x ) the vice- versa is the case. For a given e l ex 
the splitting of one of the rest two strength parameters, 
(s l +eQ l ) and (s l +e" ! ), can be decided by assuming a stan- 
dard value of symmetry energy E s (po) at normal density 
Pq. The splitting of the remaining parameter is decided 

from the value of E' s (po) = po—jf^-\p= Po - In order to de- 
cide the value of E' s (p ) we have assigned arbitrary values 
to it and calculated the EOS of NSM in each case solv- 
ing the charge neutrality and beta stability conditions. 
It is found that for a characteristic value of E' s (po) the 
asymmetric contribution to the nucleonic part of the EOS 
of NSM (that solely determines the composition of nor- 
mal neutron stars) gives stiffest behaviour that remains 
almost stationary within a small range around this value 
of E' s (po) [35] . This is referred as the universal high density 
behaviour of the asymmetric contribution of the nucleonic 
part of the EOS in NSM. We have considered this char- 
acteristic value of E' s (po) that corresponds to the stiffest 
behaviour. For the standard value of E s (p ) — 30 MeV, 
the value of E' s (po) is obtained to be 21.51 MeV for the 
EOS having 7 = 1/2 and e\ x = {e l ex + ef x ) / 3 . E' s {p ) varies 
from 20.93 MeV to 22.08 MeV as e l ex changes from to 
(e l ex + s^ x )/2 showing a small variation. Similarly E' s (p ) 
also shows a slow variation on the choice of 7 where it 
varies from 20.98 MeV to 21.70 MeV as 7 varies from 1/3 
to 2/3 (corresponding to nuclear matter incompressibil- 
ity in the range 220 to 253 MeV) for e l ex =(s l ex + sf x )/3. 
The slope parameter L — 3E' s (p n ) thus predicted for the 
EOSs considered lye within the range obtained from si- 
multaneous analysis of neutron skin thickness results in 
nuclei and neutron star properties 38,39 40 . Now, with 



the knowledge of all these nine parameters we are still left 
with two interaction parameters free for the calculation of 
finite nucleus. Here we considered to and xq of our inter- 
action in Eq.(3) as the free parameters. We determine the 
parameter to by using a simultaneous minimization along 
with the WS density distribution parameters to fit to the 
binding energy of 40 Ca nucleus. 



2.3 Determination of the parameter to and 
Wood-Saxon density distribution of nucleus 



pnucl 



The total energy of a nucleus is given by E — E" 
E Coul +E CM , where E nucl , E Coul and E CM are the con- 
tributions from nuclear, Coulomb and center-of-mass cor- 
rection, respectively. The nuclear part of the energy for an 
effective interaction is given by, 



EmCl = &i j^{r)+r n {r)]d\ 
E [?/ I Ps{r) P Ar>T\\r-r>\)d a rd*r> 



+ 2j J Ps(ry)pAryK s x '(\r-r'\)d 3 rd 3 r'], (4) 

where, the first term is the kinetic energy, second and third 
terms are direct and exchange contributions of the nuclear 
interaction; T n (r), p s (r) and p s (r,r') with s = n,p are the 
respective kinetic energy densities, densities and density 
matrices which are expressed in terms of single particle 
wave functions as, 



r(r) = 27^ 1 V^(r).V0 i (r) 
p(r) = Etx<t>*{r)Hr) 
p(r,r') = Sti4>t(r)Mr') 



(5) 



In these expressions (f>i(r) are single particle wave func- 
tions, where the subscript i denotes all the quantum num- 
bers. Instead of going into the Hartree-Fock calculation 
of single particle states we have adopted a theoretically 
transparent and numerically simplified approach where 
the density matrix expansion (DME) of Negele and Vau- 
therin [5T] is used for the density matrices in the exchange 
interaction term of the energy expression in Eq.(4). The 
density matrix under the DME can be expressed as, 



3ji(k f (R)t) 
k f (R)t 



p(R) 



35j 3 (k f (R)t).l_ 2 



2k 3 f (R) 



l^p{R) 



r{R) 



+ -k}(R)p(R)} + 



(6) 



where, j\ and j'3 are spherical Bessel functions of order 1 
and 3 respectively, kf(R) is the Fermi momentum corre- 
sponding to density p(R) at the center-of mass R of the 
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two interacting nucleons. Now choosing the Fermi momen- 
tum in the form 1421. 



Qf(R) 



5[r(fl) - \V 2 p(R)} 
3p(R) 



(7) 



reduces the DME in eq.(6) to the well known Slater ap- 
proximation p (R+i,R>-i)= 3 JlM^l p (R) of the ex- 
change term but with a modified Fermi momentum that 
accounts for the surface corrections up to second order in 
the Thomas-Fermi model. We can now express the nuclear 
part of the energy in eq.(4) as 

E nucl = J H(R)d 3 R, (8) 
where, H(R) is the energy density given by, 



H(R) = —[T n (R) + r n (R)} 



£ 2Y Ps{r) I PAr')vf(\r-r'\)d 

Zh{q a {R)t) zh(q s >{R)t) 



+Ps{r)p S '{r') 



q s (R)t q a >(R)t 



vZ (t)dH [9) 



with t being the relative coordinate and q s , s = n,p is 
the corresponding modified Fermi momentum that can be 
defined from eq.(7). The calculation of energy of a nucleus 
now requires the knowledge of density and kinetic energy 
density. Wood-Saxon density distribution, 



p{r) = 



Po 



(10) 



exp[(r — c)/a] 

is taken for simplicity that gives good description in the 
intermediate and heavy mass region. Further we make the 
approximation that p n and p p are proportional to neutron 
and proton numbers N and Z respectively. The kinetic 
energy density t s (R), s — n,p is taken to be the Thomas- 
Fermi one along with the second order correction, 



r s (R) = -k 2 sP (R) 



i \yps(R)] 

36 Ps (R) 



; V 2 p s (R), (11) 



where, k n ^{R) — [37r 2 p n ^{R)]^ is the neutron (proton) 
Fermi momentum at density p n ^(R). The WS parame- 
ters po; c and a are determined by minimizing the total 
energy, including Coulomb and center-of-correction, with 
respect to these parameters. The Coulomb energy of the 
nucleus, both direct and exchange parts, has been calcu- 
lated for the WS charge distribution. In minimizing the 
total energy we have taken a — 0.47 fm and used the nor- 
malization A = J p(r)d 3 r to express po in terms of c that 
reduces the minimization of the total energy with respect 
to c only. Now, varying to the minimization with respect 
to c is done for the experimental value of the total en- 
ergy of 40 Ca. The parameter t = 481.86 MeV fm 3 thus 
obtained for the EOS having 7 = 1/2 and (e l ex + ef x )/3 
predicts the charge radius of 40 Ca to be 3.49 fm as well 



as the binding energies and charge radii of the closed shell 
nuclei over the periodic table to an satisfactory extent as 
given in Table- 1. On varying xq from -1 to +1 the binding 
energies of N ^ Z closed shell nuclei show a slow varia- 
tion from a relatively smaller value to a higher value of 
the binding energies as compared to the results for xq = 
given in table- 1. In view of the small variation we have 
taken xq = in rest of our calculations that makes the to 
part of the interaction spin independent. With the value 
of the to thus obtained for the EOS, the binding energies 
of the proton radioactive nuclei have been calculated us- 
ing the same minimization procedure with respect to the 
WS parameters and it has been found that the results are 
reproduced within 1/2% of the experimental values for 
these nuclei. The WS density distributions thus obtained 
for these proton radioactive nuclei are used in the evalua- 
tion of the p-N interaction potentials in order to calculate 
the half-lives. 



Table 1. Binding energies B and charge radii r c of the closed 
shell nuclei. The experimental values are given besides the cal- 
culated ones inside the parenthesis. 



Nucleus 


B [MeV] 


r c [fm] 


48 Ca 


415.403 (415.991) 


3.681 (3.484) 


yu Zr 


780.928 (783.893) 


4.373 (4.272) 


x ^Sn 


1091.70 (1102.86) 


4.936 ( ) 




1628.246 (1636.446) 


5.664 (5.505) 



2.4 Nuclear part of p-N interaction potential with the 
YENI 

The p-N nuclear interaction potential given in Eq.(l) for 
the YENI in Eq.(3) becomes, 



V N (r) = VjT°(r) + V^(r) + V^ e 7» + V^(r), 

(12) 

where V^f ro contains both direct and exchange contribu- 
tions from the zero range parts of the interaction, V^^itlex) 
denotes the contribution from the finite range direct (ex- 
change) part of the interaction and V^ earr is the rear- 
rangement contribution. These various contributions are 
given by 



finite , 



V§ er °(r) = |[(1 - x Q ) Pp (r) + (2 + x )p n {r)] 
+ |f [(1 - x 3 ) Pp (r) + (2 + ^(rMj^r 



V, 



finite 



N.di 



(r) 



4ir(W+B/2-H-M/2) r e~^ r 
/i 2 L r 



f Q r' p p {r') Sl nh^r')dr' + s ™*lBZl £» r ' p p (r')e-^' dr'} 

+ W+B/2) jr r , pn{r r )smh ^ r , )dr , 



finit* 



AT,- 



2n{M-W/2+H/2-B) 
fir 
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2ir(M+H/2) 



x Jo°° ^dr' jj^j p n (R)^j^j (k(R)t)e-*dt 

V rearr {r) = _ Xs) E^M + ( 2 + X 3 )p„(r) Pp (r)] 



VS ir {r)=4ne 2 [- r' 2 Pp (r')dr' + r'p p (r')dr'} (15) 



and 



[l+fcp(r)]T+ 1 



(13) 



V^(r) = -e a (-)»4(r , )» 

7T 



(16) 



In obtaining v£™' e in Eq.(13), 
the density matrices pi(r,r'), i = 
Slater terms, 



we have approximated 
n,p by their respective 



Pi(r,r') 



3jMR)t) 
qi {R)t 



Pi(R) 



(14) 



respectively. The total Coulomb potential Vc(r) — VJj(r)- 
Vq x (r) . The Coulomb potential between the emitted pro- 
ton and the daughter nucleus is calculated from Eqs.(15) 
and (16) using the WS proton distribution of the daughter 
nucleus. 



with modified Fermi momentum q n r p ) (R) = [37r 2 p rl (p)(i?)] 1 / 3 
that can be defined from Eq.(7). The zeroth order Bessel 
function j (k(R)t) appearing in the expression of V^™* e (R) 
is a function of the wave number k(R) of the emitted pro- 
ton that contains the potential Vn(R) itself as can be 
seen from Eq.(2) and hence required to be evaluated self 
consistently. 

The p-N nuclear part of the potential, Vjy(r), in the 
cases of different proton radioactive nuclei are calculated 
from Eqs.(12,13) for a given EOS with the WS density 
distributions of the nuclei obtained from the minimiza- 
tion procedure discussed in the last sub-section. We have 
considered altogether five EOSs, two cases of different p-n 
effective mass splittings and three cases of nuclear matter 
incompressibility. In case of each of the five EOSs the pa- 
rameter to is obtained as discussed in the last subsection 
and the binding energies of the nuclei are verified to be 
reproduced within the same accuracy in each case. The 
two cases of effective mass splitting in nuclear matter at 
normal density po corresponding to the values of e l ex — 
(e' ex + e£)/6 and 
of asymmetry /3 = 



^ ex ' 

Pn-Pp 
P-n+P V 



£ex)/ 2 are shown as a function 
in Figure- 1. The results of the 



calculations of Vjy(r) for different radioactive nuclei for 
these two cases of e l ex having a given 7 value show little 
difference. In view of this insensitivity of the p-N interac- 
tion potential to the n-p effective mass splitting we have 



considered a representative value, e 



I 

ex 



)/3, in 



our subsequent calculations of p-N potentials for the three 
different cases of 7, namely, 7 = 1/3, 1/2 and 2/3, corre- 
sponding to the values of nuclear matter incompressibility 
220, 240 and 253 MeV, respectively. The results of V N (r) 
for these three EOSs are shown in Figure-2 for the case 
of 113 Cs. The difference in the results in these three cases 
are small having the characteristic behaviour of small ex- 
tension of the tail in case of lower incompressibility. We 
shall examine the effect of these variations in Vn(t) on 
the proton half-lives. 

2.5 Coulomb part of p-N interaction potential 

The direct and exchange parts of the Coulomb interac- 
tion potential of a proton with a nucleus having charge 
distribution p p (r) are given by 



3 Proton radioactivity 

In the present work, the tunneling probability of the pro- 
tons is calculated in the WKB framework. The WKB 
method has been found to be quite satisfactory for the 
a decay half life calculations and somewhat better than 
the S-matrix method [43 a . The barrier penetrability P in 
the improved WKB |44j framework for any continuous 
(rounded) potential barrier is given by 



P =!/[! + exp(K)} 



(17) 



and the decay constant by A = vPS p where S p is the spec- 
troscopic factor and the assault frequency v is calculated 
from E v = ^hv : the zero point vibration energy. The half 
life is obtained from Tx/2 — In 2/ A. The decay half life T111 
of the parent nucleus (A, Z) into a proton and a daughter 
(Ad,Zd) can, therefore, given by 

T 1/2 = [{h\n2)/(2S p E v )][l + exp{K)] (18) 

where the action integral K within the improved WKB 
approximation is given by 



K = {2 /K) / [2/i(V(r) — E v — Q)} 1/2 dr 



(19) 



with R a and Rb being its 2 nd and 3 rd turning points de- 
termined from the equations 



E{R a ) = Q + E V = E{R b ) 



(20) 



whose solutions provide three turning points. The proton 
oscillates between the first and the second turning points 
and tunnels through the barrier at R a and Rb- The zero 
point vibration energy E v is assumed to be proportional 
to Q value of the spontaneous emission of protons. For 
the present calculations, the zero point vibration energies 
used here are the same as given by Eq.(5) of Ref. [5S] 
but extended to protons and the experimental Q values 
[I] are used. The spectroscopic factor appearing in the 
denominator in Eq.(18) contribute a term —logS p to log 
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Table 2. Comparison between the measured and theoretically calculated half lives of proton emitters. The experimental Q 
values, half lives and I values are from Ref. pQ. The results of the present calculations using the YENI folded potentials are 
compared with the experimental values along with the results of DDM3Y [7] and GLDM [12| . The turning points R2=R a and 
R3=Rb are for YENI folded potentials for the case of 7=1/2 and s ex — (s ex + e"^)/3. Experimental errors in Q values Q] and 
corresponding errors in calculated half lives are inside parentheses. Asterisk symbol in the parent nucleus denotes isomeric state. 



Parent 


I 


Q ex 


R2 — Ra 


i?3 — Rb 


Measured 


YENI 


S 


expt 
p 


Op 


GLDM 


DDM3Y 


A Z 


h 


MeV 


[fm] 


[fm] 


log 10 T(s) 


log w T(s) 








logioT(s) 


log w T(s) 


iub S6 


2 


0.491(15) 


6.61 


134.30 




2.01(46) 





914 


0.999 
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4 Results and conclusion 

The half-lives in the cases of proton emitting nuclei away 
from proton drip line are calculated using the WKB bar- 
rier penetration method. The experimental Q values to- 
gether with their uncertainties are considered in calculat- 
ing the penetration probabilities in different proton emit- 
ting nuclei. The nuclear part of the p-N interaction poten- 
tial is calculated using the semiclassical approximation up 
to second order for the kinetic energy densities as well as 
for the density matrices. The DME used for the density 
matrices along with the modified Fermi momenta takes 
care of the surface corrections up to second order. The 
WS density distributions of the nuclei predicted by the 
interaction are used to calculate the nuclear and Coulomb 
parts of the p-N interaction potential. The direct part of 
the nuclear potential is evaluated exactly whereas the ex- 
change part is approximate up to the second order correc- 
tion of the density matrix expansion. The exchange part 
of the nuclear potential is evaluated self consistently. The 
interaction potential thus obtained for the YENI effec- 



tive interaction in case of each nucleus is used to cac- 
ulate the penetration probability. The half life is calcu- 
lated from Eq.(15) under the consideration that the spec- 
troscopic factor S p — 1. The results of different proton 
emitting nuclei are given in Table 2 for the EOS corre- 
sponding to 7 = 1/2 and s l ex = (e l ex + e^)/3 along with 
the respective Q and I values. The results of other calcu- 
lations using WKB barrier penetration are also listed in 
the same table for comparison together with the exper- 
imentally measured results. The agreement between the 
results of the present calculation and those of DDM3Y 
[IS] are good being close to the experimental values com- 
pared to the JLM [5] model. In the JLM model the WKB 
penetration probabilities are calculated from the interac- 
tion potentials by folding the JLM effective interaction 
with the densities of the nuclei obtained from the rela- 
tivistic mean field model (RMF). For the cases of 147 Tm, 

150 Lu> 156 Ta ^ 156 Ta * j 177 Tl * &nd 185^ The agreement 

of the calculated values in the present case as well as that 
of DDM3Y with the experimental results do not match 
well. These large deviations, particularly 177 Tl* and 185 Bi 
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Fig. 1. Neutron and proton effective masses [M*(k = 
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Fig. 2. The p-N potential V(r) and its nuclear part Viv(r) are 
shown as a function of distance r in case of the nucleus 113 Cs 
for DDM3Y and for three different cases of 7 = 1/3, 1/2 and 
2/3. All three results corresponds to same e ex = (e ex +e"j.)/6. 
For details see the text. 



could be brought down to reasonably close range of the 
experimental values in the GLDM model [12] calculation 
by including the spectroscopic factors calculated from the 
RMF+BCS theory. The spectroscopic factor is found to be 
greatly affected by the proton shell structure and in turn 
contains shell effect to a large extent. The uncertianty 
in the present calculation of half lives attributed to the 
spectroscopic factor can be obtained from the relation [5] , 

gexpt _ ^ wnere t t ex > i s the calculated (measured) 

y 1/2 ' 

haf life. The experimental spectroscopic factors Sp Xpt ob- 
tained from this relation in the present calculation is com- 
pared with the theoretical spectroscopic factors S* h cal- 
culated using RMF+BCS model. The agreement qualita- 
tively reproduces the general trend. It is worthwhile to 
mention here that the Sp Xpt values somewhat large com- 
pared to unity such as for the cases of 150 Lu and 156 Ta*, 
can be brought down from 2.884 and 3.724 to 2.317 and 
2.350, respectively, if instead of mean values extrema val- 
ues of measured and theoretical half lives are used. The 
discrepancies may be attributed to the fact that in the 
present calculations the shell and deformation effects are 
not considered rigorously. However, in the cases where the 
shell effects are not strong the calculated half lives matches 
with the experimental ones to a reasonable extent. 

The effect of n-p effective mass splitting on the decay 
half lives has been examined by claculationg the half lives 



for the two cases of e 



ex 



-.ul 



)/3 and (4 3 



)/2. 



In Fig.l, the neutron and proton effective masses [M*(k 
kf n!P ,po,a)/M] n ,p as a function of isospin asymmetry f3 
for the two cases of splittings of exchange strength pa- 
rameter into like and unlike channels. The results for half 
lives are almost same as expected from the results of the 
interaction potentials in these two cases those differ within 
the line width. Thus the n-p effective mass splitting in fi- 
nite nuclei on the proton decay has little effect and the 
reason may be attributed to the fact that the asymmetry 
as well as the Fermi momenta involved are small. On the 
other hand, the effect of the variation of nuclear matter 
incompressibility show observable effect on the decay half 
lives. With decrease in the value of K(po), the decay half 
lives decreases. By decreasing K(po) from 240 MeV (cor- 
responding to 7 = 1/2) to 220 MeV (7 = 1/3) the results 
of the calculated decay half lives decrease on the average 
byl0%. In case of 113 Cs, the change in logioT(s) is from - 
5.55 (3.12 psec) to 5.51 (2.83 psec) as 7 decreases from 1/2 
to 1/3. Similarly as K(po) is increased from 240 MeV to 
253 MeV by increasing 7 from 1/2 to 2/3, the calculated 
half lives increase on the average by 7.5 % and in 113 Cs 
logioT(s) is increased to -5.47 (3.39 psec). In Fig. 2, the 
nuclear interaction potentials V^(r) for these three EOSs 
of YENI have small differences in the tail region where 
the second turning point is located. Relatively slower rate 
of vanishing of the attractive nuclear potential V/v(r) is 
observed in the tail region for the EOS corresponding to 
a lower value of K(po) resulting in the shift of the po- 
tential barrier to higher distance. Accordingly the second 
turning point in this case will shift to a relatively higher 
distance compared to EOS corresponding to higher value 
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of K(po). In case of 113 Cs in Fig. 2 the position of the sec- 
ond turning point has decreased from 6.84 fm to 6.72 fm 
as K(p ) has increased from 220 MeV to 253 MeV. This 
shift in the position of the second turning point is solely 
determined by the nuclear part of the interaction potential 
as the Coulomb and the centrifugal parts are same for all 
these EOSs. The position of the third turning point in all 
the cases is solely determined from the Coulomb interac- 
tion potential that gives the same result for the different 
EOSs considered. Thus the width of the potential barrier 
in case of EOS corresponding to a lower value of incom- 
pressibility decreases in comparison to the EOS having 
higher incompressibility resulting in the higher penetra- 
tion probability. The values of our potential and that of 
DDM3Y [46 in the tail region in Fig. 2 are close (more 
so when Coulomb and centrifugal potentials are added) 
and hence the agreement in the results for the half lives 
as given in Table. 2. It shows that the predictions of the 
proton decay half lives in different models using the WKB 
penetration crucially depends on the nuclear potential in 
the tail region that determines the position of the second 
turning point and hence the penetration probability. 
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